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Abstract 

Inequalities for QCD functional integrals are used to establish that up 
to certain technical assumptions, the high temperature chirally restored 
phase of QCD is effectively symmetric under U(Nf) x U(Nf) rather than 
SU(Nf) x SU(Nf). If these assumptions are correct, there are no effects due 
to anomalous breaking of U(1)a on correlation functions in this phase. 
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One of the most important features of QCD is that it has an approximate SU(Nf) x 
SU(Nf) chiral symmetry which is spontaneously broken. For the purposes of this letter, it 
will be assumed that the symmetry is exact or, more precisely, that corrections due to finite 
quark masses are small and can be handled via chiral perturbation theory. If one were to 
study QCD above some critical temperature this symmetry will be restored. The nature of 
this restored phase is of more than purely theoretical interest since ultrarelativistic heavy 
ion collisions are expected to lead to thermalized regions of space with a temperature above 
T c . For simplicity in the present discussion it will be assumed that, Nf = 2. 

This letter addresses the question of what can be learned about the chirally restored phase 
directly from QCD via purely analytic means. If one makes certain technical assumptions, 
one can deduce nontrivial — and rather surprising — things about the nature of this phase 
by exploiting QCD inequality techniques similar to those used by Wiengarten, Vaffa and 
Witten in studies of QCD at T = 0. 

In particular it will be shown that if a certain set of zero measure does not afflict the 
functional integral then, above T c , the phase is effectively symmetric under {7(2) x U (2) in 
the sense that operators can be classified into multiplets associated with representations of 
{7(2) x i7(2) and that correlation functions of operators in a given multiplet are identical. 
This means, for example, that the two-point correlation function in the ir channel is degen- 
erate with the correlation function in the rj channel. This is surprising since the U (2) x {7(2) 
symmetry of the QCD lagrangian is broken by the U(1)a anomaly to SU(2) x SU(2). More- 
over, the U(1)a anomaly is at the operator level and thus the anomaly exists independent 
of temperature. As stressed by 't Hooft 0, the anomaly may be thought of as providing a 
mechanism for explicit (as opposed to spontaneous) symmetry breaking. Thus, it seems a 
priori implausible that restoration of chiral SU(2) x SU(2) should imply invariance under 
{7(2) x {7(2). Indeed, in a classic early review of the subject of instantons and the U(1)a 
problem, Coleman |J asserts precisely the viewpoint that the U(1)a symmetry remains 
broken above the restoration temperature with the symmetry breaking decreasing at high 
temperatures as a power of 1/T. Moreover Meggiolaro Q has recently constructed a model 
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motivated by lattice calculations of the topological susceptability H in which U{1)a remains 
broken above the chiral restoration temperature. 

The idea that the chirally restored phase is symmetric under U(2) x U(2) rather than 
SU{2) x SU{2) is not new. Shuryak raised this possibility previously [/J. The present 
approach is novel, however, in that it derives this result formally from the QCD functional 
integral. Before discussing the present derivation, it is useful to review Shuryak's arguments. 
One argument is based on lattice calculations of screening masses which purport to show 
that above T c , the ir and a screening masses are degenerate (within numerical noise) |§. The 
calculated "<r" correlation functions only included the quark-line connected part. However, 
this quark-line connected part is the entire correlator in the scalar-isovector (S) channel. 
Thus, the lattice calculations indicate a degeneracy between the n and S screening masses. 
The it and the 5 do not belong to the same SU{2) x SU{2) multiplet; they do, however, 
belong to the same U(2) x £7(2) multiplet. 

Shuryak also argues for U(2) x U (2) restoration from the instanton liquid model P,|10 



Recall that the solution of the U (1)^ problem requires both the U(1)a anomaly and the con- 
tribution of nontrivial topological configurations [[J ; topology is necessary for the anomalous 
violation of U(1)a symmetry to have physical manifestations. In the instanton liquid model 
PJlO[| instantons provide the only source for these configurations. As has long been known, 



a finite density of instantons leads to chiral symmetry breaking [LlJ. In the instanton liquid 
model, instantons also provide the only source of SU(2) x SU(2) chiral symmetry breaking. 
Thus, in this model, the same mechanism is responsible both for SU(2) x SU(2) chiral sym- 
metry breaking and for allowing the U(1)a symmetry breaking due to the anomaly to have 
physical consequences. In such a model, if one were in a phase in which SU(2) x SU(2) chi- 
ral symmetry were unbroken, it would follow that instanton effects are be turned off. This, 
in turn, suggests that the anomaly will not have physical effects in any of the correlation 
functions: all observables will behave as though the phase is U(2) x U(2) symmetric. The 
mechanism responsible for this in the model is believed to be the condensation of instantons 
and anti-instantons into topologically neutral "molecules" [|l!| . 



The present analysis is based on properties of the QCD functional integral. The essential 
physics is best understood from the quark propagator in a given gluon background field. 
When the propagator is written in terms of a spectral representation, all U(1)a violating 
effects come from eigenmodes in the neighborhood of zero virtuality; i.e. A = modes, 
where the Dirac eigen equation is pipj = iXjipj- While it is not immediately obvious how to 
establish in general that all U(1)a violating amplitudes come from the region of A = 0, it is 
easy to establish for given U (I) a violating amplitudes by studying the spectral representaton 
of the propagator in the context of the functional integral. Moreover, it is easy to prove 
that the density of states of the Euclidean Dirac operator, Jj) at A = 0, is zero for any gauge 
configurations with boundary conditions consistent with a temperature greater than the T c 
(excluding, perhaps, a set of zero measure). Thus, one can show that these U(1)a violating 
amplitudes vanish. 

The fact that above T c all gauge configurations yield a vanishing density of states at 
zero virtuality can be seen quite transparently. The chiral condensate (qq) is related to the 
density of states at zero averaged over gluon field configurations JT3]: (qq) = —tt(pa(X = 
0)) where pA is the density of states in a given background gluon field configuration and 
( ) indicates averaging over the gluon field configurations weighted by e~ SYM Det[p — m]. 
This applies to the finite temperature case provided the average over gluons only includes 
configurations periodic in Euclidean time with a periodicity of /3 — 1/T and the fermion 
determinant is evaluated for antiperiodic configurations. Above T c , ((qq))r = 0, implying 
that ((pa(0)))t = (where the double bracket indicates a thermal average). However, 
PaW is a density; accordingly it is positive semi-definite (i.e. > 0): Pa(A) > 0. Moreover 
the weighting function e~ YM Det[p — m] is also positive semi-definite [|IJ. An averaged 
quantity which is never negative cannot have an average of zero unless the quantity is zero 
for all configurations (except, perhaps, a set of measure zero): Pa(0) = for configurations 
consistent with the boundary conditions for T > T c . 

Before discussing how this works out in specific cases, it is worth stressing the generality of 
the result. It depends only on the fact that Pa(0) goes to zero above the phase transition for 
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all configurations and that U(1)a violating amplitudes come from modes in the neighborhood 
of A = 0. It does not depend on the detailed mechanism which generates a nonzero Pa(0) 
below T c . 

To make the discussion concrete, consider the two-point correlation function of scalar 
and pseudoscalar quark bilinears. There are four distinct operators: pseudoscalar-isovector, 
iqj^rq, (the n channel); scalar- isoscalar, qq (cr); pseudoscalar-isoscalar, iq^^q, (r]'); and 
scalar-isovector, qrq (5). These bilinears are denoted as J n , J a , J v > and J$- The o and it 
form a distinct SU{2) x SU(2) multiplet from the S and rj'. Under U(2) x U(2) however 
they are all part of a single multiplet. 

The thermal two-point correlation function Il(x) of two equal-time quark bilinear oper- 
ators at fixed temperature, T, is defined by 

Hj(x) = «J(x)J(0)» T - «J(x)» T «J(0)» T (1) 

where the double braces indicates thermal average and J(x) = q(x)Tq(x) and T is a matrix 
in Dirac and flavor space. One can write this as a Euclidean functional integral: 

n 7 (x) = ~ f D[A] e- s ™Det[p-m q ] (2) 

Z JT 

x tr[£4x,0)r<S A (x,0)r] -tr[S A (x,x)r]tr[,S A (0,0)r] 

where the subscript T indicates the finite T boundary conditions (periodic in A); Z is the 
partition function; Sym is the Euclidean action of the Yang-Mills field; the Det indicates a 
functional determinant with the fermion modes satisfying antiperiodic boundary conditions; 
SU(x, y) is the Euclidean space quark propagator in the presence of a background gauge 
field A; and the traces are over color, flavor and Dirac spaces. 

A finite quark mass m q is included in the previous expression — it will be sent to zero 
only at the end of the calculation. For technical reasons it is simpler to work in a box of 
finite volume V (which makes all of the modes discrete) and to let the volume of the box go 
to infinity at the end of the problem. The ordering of these two limits is critical. One must 



take the V — > oo limit before taking the chiral limit. [14 



There are two distinct contributions to this functional integral: a term with a single trace 
and a term with two traces. They are the quark-line connected and quark-line disconnected 
pieces, respectively. If the up and down quark masses are equal (as is assumed here), 
[«Sa(x, 0), r] = 0, and the connected piece of an isoscalar correlator (e.g. the a channel) is 
identical to the connected piece of an isovector correlator with the same spatial quantum 
numbers (e.g. the 5 channel). 

The difference between the a and 5 correlation functions, U a (x) — Hs(x), is U(1)a vi- 
olating. As noted above, in the functional integral this difference comes entirely from the 
quark-line disconnected piece, 

n CT (x) - n 5 (x) = \\ t D[A}e' s ^Det[p-m q }tr[S A (^^)]tT[S A (0,0)} . (3) 

If it can be shown that above the SU (2) x SU (2) chiral restoration temperature 

tr[^(x,x)]~OK) (4) 

for all gauge configurations consistent with the boundary conditions, then it follows that 
tr[5 < A(x, x)]tr[S' J 4(0, 0)] ~ 0(m 2 q ) for all gauge configurations and thus the weighted average 
over gauge configurations will also be 0(m 2 q ) from which eq. (|3|) implies 

n CT (x) - ILi(x) ~ 0(m\) . (5) 

This in turn implies that in the chiral limit of m q — > II (7 (x) — 11,5 (x) ~ * 0. That is, this 
U(1)a violating matrix element vanishes. 

If the validity of eq. (|j) is established, then one has proven that this U(1)a violating 
amplitude vanishes. To begin use a spectral representation for Sa' 

q i x ^ ^(x)^j(y) 

where the modes are eigenmodes of the Dirac operator. ^From the fact that {75, p} = 0, 
it follows that if ipj is an eigenmode with eigenvalue then 75^ is an eigenmode with 
eigenvalue — This in turn implies 



tr[5 A (x,x)]=^ jf— -j • (7) 

It is apparent from eq (|?P that, as advertised, in the limit of m q — > 0, contributions to 
tr[S^(x, x)] come entirely from modes near A = 0. More significantly, given the standard 
convention that the quark mass is positive, then trfS^x, x)] < for any gauge configuration. 

Above T c , we know that chiral condensate, ((qq))T, vanishes, or to be more precise is 
order m q and vanishes when the chiral limit is taken. The chiral condensate can be written 
clS db functional integral: 

N f {{qq(x))) T =^J t D[A] e" 5 ™ Det[p- m q \ tr[S A (x, x)] ~ -0(m q ) (8) 

At this stage, it is worth recalling that e~ SYM Det[p— m q ] is positive semi-definite for all gauge 
configurations while tr[SU(x, x)] is negative semi-definite. Thus, the integrand in eq. (§) is 
negative semi-definite. This means that there can be no cancellations in the integral — the 
only way that the integral can be 0(m q ) is if the contributions from all gauge configurations 
are 0(m q ) (except perhaps from a fraction of configurations which goes to zero in the chiral 
limit). Thus eq. (|J) has been shown to be true for all gauge configurations contributing to 
the functional integral except for contributions which become a set of measure zero in the 
chiral limit. Assuming this set of measure zero can be safely ignored in the evaluation of 
eq. ([D — an issue which will be discussed at the end of this letter — one concludes that since 
eq. (JD is true so is eq. ([5]); thus in the chiral limit of m q — > these a and 5 correlators are 
identical. 

Having established this, it is immediately obvious that U n (x), U v '(x), IL^x), and ILj (x) 
must all be identical above T c in the m q — > limit. SU(2) x SU(2) chiral restoration implies 
that n^x) = n CT (x) and ITy (x) = ITj (x), while eq. (|^) implies that rbj(x) = n CT (x) - 
all members of the U(2) x U(2) multiplet are identical. Although from this argument it is 
clear that ITy(x) = Il vr (x), it is useful to demonstrate this directly from functional integral 
inequalities as it demonstrates a technique which is useful for studying other multiplets. 

The functional integral for this difference can be written as 
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^(x) - ILy(x) = -J J D[A]e- 5YM Det[|)-m (? ]tr[5 A (x,x)7 5 ]tr[^(0,0)75] • (9) 

The first step in proving that II T (x) — Il f? /(x) goes to zero above T c is to show that 

|tr[S A (x,x) 75 ]| < |tr[^(x,x)]| (10) 

for any gauge configurations. This is easily established using the spectral decomposition 
of the propagator and the fact that ■0j(x)(l + 75) 2 ^j(x) > for any ipj. By comparing 
eq. (|D with eq. (|) and using eq. ( |10|) and the fact that e~ SYM Det[p— m q ] is positive semi- 
definite, one sees that |n„.(x) — II^(x)| < (n^x) — (x) | . Since the right-hand side of 
this inequality goes to zero, the left-hand side does as well, and thus the degeneracy of the 7r 
and rj' channels above T c has been demonstrated directly from the QCD functional integrals. 

The technique used to establish that ir and rf correlation functions are identical above 
the phase transition by showing that the absolute value of their difference is less than or 
equal to jll^ — Hg\, can be immediately generalized for other channels. In this way, one 
can show that vector and pseudo-vector, isovector and isoscalar (i.e. the ui, p, f\ and a{) 
correlation functions are all identical above T c . Again this is an identification of U (2) x U(2) 
symmetry since only the p and a\ are connected by SU (2) x SU(2) chiral symmetry. The 
same method allows one to show that the tensor and pseudo-tensor, isoscalar and isovector 
correlators are identical above T c . 

There is a loophole in the demonstration of the U(2)xU (2) nature of the chirally restored 
phase given above. In particular, it was assumed that contributions to the functional integral 
in eq. (||), which were a set of measure zero in the m q — > limit, do not contribute to the 
functional integral in eq. (0). By inspection, it is clear that so long as tr[S < A(x, x)] is finite 
for all gauge configurations (after a gauge invariant and SU(2) x SU(2) chiral invariant 
ultraviolet regularization), then the set of measure zero cannot affect the functional integral 
in eq. (|3|). To discuss a set of measure zero with infinite contributions it is sensible to first 
introduce an infrared cutoff regulator, e where e — > corresponds to the infinite volume and 
m q — > limits with Vm? q — > oo. The loophole in the general argument is that there could be 
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configurations for which tr[S^(x, x)] ~ e -1 / 2 which have a weight proportional to e. In such 
a case, ((qq(x))) T ~ e 1 ^ 2 which vanishes in the e — > limit while {n CT (x) — ^(x)} ~ 0(1). 
Thus, there is apparently the possibility that the chiral condensate vanishes as the regulator 
goes to zero while the U(1)a violating amplitude, {n CT (x) — ILj(x)}, does not. 

In summary, up to the loophole discussed above, it has been shown directly from the 
QCD function integral that above T c , the correlation functions for quark bilinears in a 
given (7(2) x (7 (2) multiplet are identical. This indicates that the phase is invariant under 
U(2) x U(2) rather than SU(2) x SU(2). The anomalous U(1)a breaking does not split the 
£7(2) x U{2) multiplets because the effects of the anomaly occur entirely through the quark- 
line disconnected parts of correlation functions and, in the m q — > limit, the quark-line 
disconnected parts contribute only due to the the modes near A = 0. Above T c the density 
of states at A = goes to zero and the anomaly ceases to play a role. 
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